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Scope  and  Method  of  Study:  The  principal  goal  of  the  study  was  to  evalu¬ 
ate  an  analytical  procedure  for  predicting  progressive  collapse  in 
damaged  complex  structures.  A  structure  was  modeled  for  analysis 
by  the  finite  element  method  using  relatively  large,  simple  ele¬ 
ments.  There  was  little  or  no  refinement  of  mesh  size  in  areas  of 
initial  damage  or  damage  propagation.  A  method  was  developed  for 
determining  and  applying  allowable  stresses  to  help  compensate  for 
the  absence  of  model  detail.  Stress  results  of  a  finite  element 
analysis  were  examined  by  a  computer  post-processor  program  written 
for  this  study  to  make  selective  changes  to  the  finite  element 
model.  The  modified  model  was  analyzed  using  the  finite  element 
method  and  the  procedure  was  repeated  in  an  iterative  fashion  to 
predict  progressive  collapse.  Analytical  results  were  compared  to 
experimental  test  data  to  determine  the  validity  of  the  analytical 
procedure.  —  -  - 

Findings  and  Conclusions:  The  analytical  procedure  provided  a  relatively 
economical  method  for  predicting  progressive  collapse  in  a  complex 
structure.  Evaluation  of  a  complex  structure  subjected  to  three 
initial  damage  conditions  showed  acceptable  correlation  between  ex¬ 
perimental  and  analytical  results.  The  method  of  determining  appro¬ 
priate  allowable  stresses  was  general  enough  to  apply  to  a  wide 
range  of  materials  and  structures.  The  procedure  proved  to  be  an 
economical  estimating  tool  for  predicting  residual  structural  strength 
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CHAPTER  I 


INTRODUCTION 

Many  structures  are  susceptible  to  progressive  collapse,  a  chain  re¬ 
action  type  of  failure  following  damage  to  a  relatively  small  portion  of 
the  structure.  The  more  specialized  a  structure  is,  the  more  vulnerable 
it  is  to  progressive  failure  largely  because  it  is  designed  to  resist 
fewer  possible  loading  conditions.  As  efforts  increase  to  optimize  de¬ 
signs  within  acceptable  factors  of  safety,  the  risks  of  initiating  pro¬ 
gressive  collapse  through  relatively  minor  localized  damage  also  increase. 
An  ability  to  predict  analytically  the  response  of  a  damaged  structure 
would  therefore  be  beneficial. 

Although  progressive  collapse  is  normally  associated  with  high-rise 
buildings,  interest  in  it  is  not  limited  to  conventional  civil  engineer¬ 
ing  applications.  The  Department  of  Defense  needs  the  capability  to  pre¬ 
dict  the  residual  strength  of  battle-damaged  aircraft  and  to  know  the 
role  of  progressive  collapse  in  that  setting.  Specifically,  the  Depart¬ 
ment  of  Defense  Joint  Technical  Coordinating  Group  for  Munitions  Effec¬ 
tiveness  is  interested  in  the  post-damaged  capabilities  of  potentially 
host i le  ai rcraft. 

In  pursuit  of  its  interest,  the  Group  provided  research  funds  and 
three  F-84F  aircraft  wings  for  this  study.  The  goal  was  to  evaluate  a 
potentially  versatile  method  for  predicting  progressive  collapse  in  air¬ 
craft  structures.  The  method  was  to  be  verified  by  experimental  testing. 
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Desirable  characteristics  to  be  imparted  to  the  method  would  be  relative 
simplicity  in  preparing  for  its  use  and  relative  ease  and  economy  in  its 
appl i cation. 

The  finite  element  method  was  the  fundamental  tool  for  determining 
stresses  within  the  wing.  In  this  report  most  discussion  of  the  finite 
element  method  is  of  a  general  nature.  The  NASTRAN  (National  Aeronauti¬ 
cal  and  Space  Administration  Structural  Analysis)  program  was  selected  to 
apply  the  finite  element  method  because  of  its  versatility  and  its  wide¬ 
spread  availability  in  both  industry  and  the  defense  community.  The  read¬ 
er  is  assumed  to  be  familiar  with  the  finite  element  method  in  general. 
Where  reference  to  specific  program  characteristics  is  essential,  a  basic 
familiarity  with  NASTRAN  is  also  assumed. 

Finally,  the  sponsor  of  this  research  is  interested  in  the  effective¬ 
ness  of  munitions  in  destroying  combat  aircraft.  Consequently,  any  con¬ 
servative  assumption  is  one  which  tends  to  give  the  structure  more  strength 
than  actually  exists.  This  definition  of  conservative  is  used  throughout 
the  study.  Caution  must  be  exercised  in  directly  extending  the  results  of 
this  study  to  more  conventional  applications.  In  such  use  the  assumptions 
of  this  study  would  become  unconservative. 
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CHAPTER  II 


LITERATURE  REVIEW 

2.1  General 

Many  papers  have  addressed  the  topic  of  progressive  collapse  of  dam¬ 
aged  structures,  but  only  one  has  provided  a  general  quantitative  method 
of  analysis  (1).  The  following  sections  summarize  the  published  papers 
while  the  last  section  details  the  one  general  approach. 

2.2  Qualitative  Analysis 

Most  studies  of  progressive  collapse,  as  applied  to  structures  con¬ 
ventionally  associated  with  Civil  Engineering,  fit  into  three  categories. 
The  first  addressed  a  need  to  predict  statistically  the  frequency  and 
severity  of  damaging  events  such  as  vehicle  impact  or  explosion  (2  through 
8). 

Another  category  was  the  qualitative  analysis  of  a  structure's  abil¬ 
ity  to  resist  damage  or  to  develop  alternate  load  paths  around  damaqe. 
Typical  tonics  of  discussion  included  catenary  action  of  slabs,  beam  ac¬ 
tion  of  adequately  tied  cei 1 i ng-wal 1 -f loor  systems  acting  as  wide  flange 
sections,  and  the  in-plane  arching  of  walls  over  damage  (4,  7,  9  through 
H). 

A  third  category  was  an  effort  to  develop  codes  which  mate  the  first 
two  areas  into  economically  and  socially  acceptable  guidelines  for  design 
and  construction  (15  through  23).  Additionally,  a  research  workshop  was 
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conducted  in  1975  to  evaluate  present  knowledge  of  the  progressive  col¬ 
lapse  phenomenon  and  to  identify  areas  requiring  further  study  (24). 

2.3  Quantitative  Analysis  of  Building  Structures 

A  smaller,  fourth  category  addressed  the  need  to  evaluate  quantita¬ 
tively  the  behaviors  occurring  during  a  progressive  collapse.  Several 
studies  have  been  completed,  but  most  have  considered  only  two-dimensional 
problems  and  most  have  required  extensive  analyst  interactive  involvement 
(25  through  29). 

Smith  and  Epstein  (30)  developed  a  three-dimensional  method  to  ana¬ 
lyze  the  progressive  collapse  of  a  space  truss  roof.  Their  approach  used 
the  finite  element  method  to  determine  structural  member  stresses.  As  a 
member  approached  its  buckling  load,  predetermined  for  every  member  in 
the  structure,  the  member  was  replaced  by  opposite  equal  forces  repre¬ 
senting  post-buckling  strength.  The  method  did  provide  a  three-di mensi on- 
al  analysis  but  was  limited  exclusively  to  buckling  related  failures.  It 
was  inappropriate  for  structures  in  which  other  failure  modes  share  equal 
importance  or  are  dominant. 

2.4  Quantitative  Analysis  of  Aircraft  Structures 

The  military's  need  to  predict  the  behavior  of  damaged  aircraft  has 
precipitated  several  papers  of  interest.  Venkayya  (31)  outlined  an  em¬ 
pirical  iterative  procedure  in  1978  for  determining  the  residual  strength 
of  damaged  structures.  The  displacements  and  decomposed  stiffness  matrix 
of  an  undamaged  structure  were  combined  with  a  sparse  negative  stiffness 
matrix  representing  damage.  The  result  was  an  iteratively  derived  second- 
order  Taylor  series  approximation  of  the  response  of  the  damaged  structure. 
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The  method  appeared  suitable  for  economic  evaluation  of  initial  structure 
response  to  several  different  damage  conditions.  However,  when  the  me¬ 
thod  is  applied  to  progressive  collapse  analyses,  problems  surface  as 
component  failure  progresses  toward  collapse.  Solution  convergence  times 
become  unacceptably  slow  and  convergence  criteria  become  increasingly 
difficult  to  establish. 

In  1976,  Heard  (1)  proposed  for  the  Air  Force  Armaments  Testing  Lab¬ 
oratory  (AFATL)  a  method  of  structural  modeling  and  analysis  for  progres¬ 
sive  collapse  in  aircraft  structures.  That  method,  referred  to  in  this 
study  as  the  AFATL  method,  appeared  to  be  the  most  promising  general  ap¬ 
proach  to  a  quantitative  analysis  of  progressive  collapse.  The  next  sec¬ 
tion  presents  this  method  in  some  detail. 

2.5  Analysis  Method  Background 

The  structure  being  evaluated  must  be  represented  as  a  computer  model 
for  finite  element  analysis.  Because  the  method  requires  many  iterative 
analyses  to  trace  the  progressive  collapse  phenomenon,  economy  urges  the 
use  of  the  largest,  simplest  elements  which  still  describe  the  basic  geo¬ 
metry  of  the  structure  and  provide  adequate  precision  to  permit  a  stress- 
based  analysis.  A  principal  feature  of  the  AFATL  method  is  that  little 
or  no  refinement  of  the  model  occurs  in  the  area  of  damage.  This  feature 
aids  the  economy  of  the  method  but,  because  the  large  elements  mask 
stress  concentrations,  the  method  must  include  compensating  techniques. 
Heard  employed  two  such  techniques  which  are  described  later. 

A  load  was  applied  to  the  model  and  the  resulting  stresses  were  exam¬ 
ined  in  search  of  overstressed  elements.  An  overstressed  element  was  one 
whose  stresses  exceeded  predefined  limiting  values.  A  solitary  over- 
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stressed  element  was  removed  from  the  model  as  having  failed.  If  more 
than  one  overstressed  element  occurred  grouped  together,  only  the  most 
severely  stressed  element  of  the  group  was  removed.  This  technique 
helped  represent  crack  propagation  in  a  model  composed  of  large  elements 
and  was  supported  by  studies  of  Sih  and  Hartranft  (28). 

Reducing  the  values  of  limiting  stresses  for  elements  bordering  dam¬ 
age  was  the  second  technique  to  compensate  for  loss  of  stress  concentra¬ 
tions  around  crack  tips  at  the  edges  of  damage.  Thus  the  computed  stress 
in  an  element  bordering  damage  might  produce  element  failure  while  a 
similarly  stressed  element  away  from  damage  remained  intact.  Using  dif¬ 
ferent  values  for  limiting  stresses  complicated  the  process  of  selecting 
which  element  to  fail  in  a  group  of  overstressed  elements.  The  most 
severely  stressed  element  could  not  be  determined  through  a  direct  com¬ 
parison  of  the  magnitudes  of  element  stresses. 

After  the  failed  elements  were  removed,  the  modified  model  was  again 
analyzed  and  the  procedure  was  repeated.  Iterations  continued  until  the 
model  could  sustain  some  desired  maximum  load,  or  until  failure  occurred. 
This  latter  condition  was  sometimes  determined  subjectively  by  evaluating 
the  structure's  displaced  shape  rather  than  by  its  residual  load-carrying 
capacity. 


CHAPTER  III 

SPECIFIC  OBJECTIVES 

The  method  proposed  by  Heard  appeared  to  be  a  versatile  approach  for 
the  quantitative  analysis  of  progressive  collapse.  To  increase  the  accep¬ 
tability  of  the  method,  however,  four  areas  were  identified  as  objectives 
for  further  study. 

Validation  of  the  method  was  perhaps  the  most  important  objective. 

Due  to  an  absence  of  actual  aircraft  wings  which  his  model  represented, 
Heard  was  unable  to  substantiate  with  actual  test  data  the  value  of  his 
work.  The  first  phase  of  this  study  was  a  laboratory  test  program  which 
provided  data  for  evaluation  of  analytical  results.  Tests  of  three  F-84F 
aircraft  wings  measured  structural  performance  under  different  damage  and 
1 oad  comb i na  t i ons . 

In  the  previous  study,  only  one  combination  of  elements  was  used  for 
modeling  the  aircraft  wing  structure.  A  comparison  of  several  element 
combinations  was  made  in  search  of  the  best  selection  of  model  elements. 

The  actual  application  of  the  method  required  a  great  amount  of 
manual  data  analysis  for  each  iteration.  A  large  number  of  elements  had 
to  be  checked  and  compared  to  limiting  stress  values.  The  relative  loca¬ 
tions  of  overstressed  elements  had  to  be  determined  and  caution  applied 
to  remove  the  appropriate  element.  Finally,  removal  of  failed  elements 
required  modifications  of  the  model.  In  addition  to  removing  the  failed 
elements,  modification  included  reducing  limiting  stresses  for  elements 
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bordering  the  newly  propagated  damage.  Automating  the  application  of  the 
method  was  desirable  to  reduce  both  time  and  expense  for  a  complete  anal- 


The  final  area  to  address  was  the  appropriate  values  for  limiting 
stresses.  Heard  used  two  levels  of  limiting  stresses:  material  ultimate 
strength  for  elements  away  from  damage,  and  material  yield  strength  for 
elements  bordering  damage.  A  more  sophisticated  determination  of  limit¬ 
ing  stresses  had  the  potential  for  returning  more  realistic  results. 

These  four  areas, 

1.  Comparison  of  analytical  and  test  results 

2.  Comparison  of  modeling  elements 

3.  Automation  of  the  method 

4.  Determination  of  limiting  stress  values, 
became  the  specific  objectives  of  this  study. 


CHAPTER  IV 


EXPERIMENTAL  TEST  PROGRAM 

4.1  General 

A  major  objective  of  this  study  was  to  use  actual  test  data  as  a 
standard  for  evaluating  analytical  results.  Three  F-84F  aircraft  wings 
were  tested,  each  with  a  different  damage  and  load  combination.  This 
chapter  contains  descriptions  of  specific  damage  and  loads  and  of  the 
general  test  procedure.  Appendix  A  contains  diagrams  showing  strain  gage 
locations  for  the  various  tests. 

4.2  Specific  Test  Descriptions 

The  F-84F  aircraft  wing  is  a  two-spar,  semi -monocoque  structure. 

For  general  reference,  Figure  1  illustrates  the  upper  and  lower  wing  sur¬ 
faces  and  the  wing  structural  frame.  All  three  wings  were  mounted  upside 
down  for  testing;  however,  the  terms  "upper"  and  "lower"  refer  to  the 
wing's  upper  and  lower  surfaces,  not  to  their  physical  orientation  for 
the  tests.  For  all  tests,  the  landing  gear  and  gear  doors,  flaps,  and 
ailerons  were  removed. 

Test  1  consisted  of  severe  damage  to  the  upper  half  of  the  front 
spar  as  shown  in  Figure  2.  A  load  applied  to  the  front  spar  produced  a 
failure  with  bending  as  the  predominant  behavior.  The  damaged  area  was 


stressed  in  tension. 


(a)  Upper  Surface 


Test  2,  Figure  3,  was  to  measure  behavior  with  a  significant  amount 
of  torsion  present.  The  rear  spar  was  completely  severed  and  the  load 
was  applied  to  the  rear  spar.  The  result  was  a  combination  of  bending 
and  torsion  in  the  front  spar.  This  wing  could  not  be  failed  within  safe 
limits  of  the  laboratory  loading  apparatus.  Consequently,  three  loading 
trials  were  performed  on  this  wing  and  designated  Tests  2A,  2B,  and  2C. 
Each  test  had  slightly  modified  damage  to  the  skin  adjacent  to  the 
severed  spar.  These  were  efforts  to  initiate  tearing  of  the  skin  over 
the  wheel  well  area;  however,  no  propagation  of  that  damage  occurred. 

Test  3  was  an  attempt  to  represent  more  closely  the  damage  which 
could  occur  from  a  shaped-charge  missile  warhead.  Figure  4  shows  a  Si- 
inch  wide  strip  of  material  removed  from  the  lower  wing  surface.  All 
skin  was  removed  from  the  strip,  which  extended  from  the  rear  spar  to  the 
leading  edge.  The  lower  rear  spar  cap  was  removed  but  the  web  was  left 
intact.  The  portion  of  the  lower  front  spar  cap  extending  from  the  web 
toward  the  trailing  edge  was  also  removed.  The  load  applied  to  the  rear 
spar  put  the  damaged  surface  into  compression. 

The  residual  strength  of  this  wing  also  exceeded  the  safe  capacity 
of  the  loading  equipment.  A  variation  of  this  test,  designated  Test  3B, 
included  further  damage  to  the  front  spar  cap.  Half  the  width  of  the 
lower  spar  cap  extending  from  the  web  toward  the  leading  edge  was  remov¬ 
ed.  A  U-inch  width  of  spar  cap  remained  extending  from  the  rear  face 
of  the  web  toward  the  leading  edge.  This  additional  damaae  led  to  com¬ 
plete  structural  failure. 

4.3  Wing  Support  System 

Each  wing  soar  root  mounted  into  a  support  structure  as  illustrated 
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Figure  4.  Damage  for  Test  3 


in  Figure  5.  Pins  secured  the  wings  in  the  support  structures  in  the 
same  manner  as  the  wings  had  been  attached  to  aircraft  fuselages.  The 
support  structures  were  extremely  rigid  compared  to  the  wings  so  that  no 
appreciable  deformation  occurred  within  the  supports  themselves.  Three 
transducers  supported  each  T-shaped  support  structure,  permitting  measure 
ments  of  vertical  reaction  forces  and  reaction  moments  about  two  perpendi 
cular  horizontal  axes. 

The  wing  spar  roots  were  aligned  with  the  support  structures  and 
pinned  into  place  within  small  tolerance;  however,  some  motion  of  the 
wing  spar  roots  with  respect  to  the  supports  was  unavoidable.  For  Tests 
2  and  3,  dial  gages  measured  relative  rotation  of  each  wing  spar  root 
about  horizontal  axes  parallel  to  and  perpendicular  to  the  root  itself. 
These  data  then  formed  the  basis  for  support  conditions  in  corresponding 
finite  element  analyses.  These  support  conditions  provided  a  better  ana¬ 
lytical  representation  of  wing  deflections;  however,  support  conditions 
assuming  no  relative  rotation  were  used  for  stress  analyses. 

Load  System 

A  movable  overhead  crane  applied  a  single  point  load  in  each  test. 
The  crane  was  self-adjusting  so  the  load  was  always  applied  vertically. 

A  cable  attaching  the  crane  to  the  wing  load  point  was  equipped  with  an 
in-line  transducer  to  permit  continuous  accurate  monitoring  of  the  actual 
load  appl i ed . 

A.  5  Deflection  Measurements 

Vertical  deflections  were  measured  at  points  along  the  front  and 
rear  spars  corresponding  to  node  points  in  the  finite  element  model. 
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Figure  5.  Wing  Support  Structure 
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Steel  scales  incremented  to  0.01  inch  were  attached  to  the  lower  wing 
surface  and  measurements  were  read  through  an  engineer's  level.  Similar¬ 
ly,  scales  were  mounted  on  the  support  structures  above  each  transducer 
to  detect  any  vertical  displacements  there.  At  the  load  point  no  scale 
could  be  attached  to  trie  wing  surface  as  at  other  locations  along  the 
spars.  Instead,  a  cloth  tape  hung  down  vertically  from  the  upper  wing 
surface  to  measure  deflections  with  respect  to  the  laboratory  floor. 

4.6  Strain  Measurements 

Strain  gages  were  mounted  to  the  wing  to  detect  changes  in  load 
paths  as  components  failed  and  to  detect  load  levels  at  which  failures 
occurred.  The  different  designs  of  each  test  and  experience  from  previ¬ 
ous  tests  led  to  slightly  different  strain  gage  placement  for  each  wing. 
Appendix  A  contains  specific  locations. 

Quarter-inch  uniaxial  strain  gages  measured  outer  fiber  strains 
along  spar  and  rib  caps.  Three-gage  rectangular  rosettes  attached  to 
selected  skin  panels  measured  panel  behavior.  Similar  rosettes  measured 
shear  in  rib  and  spar  webs  in  Test  3* 

Wings  were  first  loaded  enough  to  compensate  for  self-weight,  and 
all  gages  were  zeroed.  For  Test  1,  all  transducers  and  strain  gages  fed 
into  a  single  switch  and  balance  unit  to  measure  output.  All  other  tests 
used  a  Vi  shay  Instruments  Measurements  Group  computer-controlled  data 
acquisition  and  reduction  system.  The  System  4000  included  the  software 
program  plus  a  Controller  4220  and  two  Strain  Gage  Scanners  4270.  A  Hew¬ 
lett-Packard  9825B,  upgraded  to  9825T  capabilities,  served  as  the  Execu¬ 
tive  Control  Unit  to  complete  the  system  hardware. 
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CHAPTER  V 


FINITE  ELEMENT  MODELS 
5-1  Background 

The  fundamental  modeling  philosophy  used  by  Heard  (1)  applied  also 
to  this  study.  Rod  elements  in  combination  with  shear  panels  represented 
heavy  structural  members  such  as  spars  and  ribs.  Shear  panel  or  membrane 
elements  represented  aircraft  skin.  Skin  stiffeners  were  modeled  by  rod 
elements. 

The  specific  structure  for  this  study,  the  F-84F  aircraft  wing,  was 
also  analyzed  by  Jordan  (32,  33).  In  1976,  he  performed  a  dynamic  re¬ 
sponse  and  small  static  load  bending  analysis  of  the  wing.  Although  his 
objective  differed  from  Heard's,  he  applied  the  same  fundamental  philoso¬ 
phy  to  develop  his  model  of  the  wing.  Jordan's  model  was  the  nucleus  of 
the  models  evaluated  in  this  study  and  is  illustrated  in  Figure  6.  De¬ 
tails  of  element  numbering  are  presented  in  Appendix  B. 

The  structure's  geometry  determined  the  size  of  the  elements.  Inter¬ 
sections  of  spars  and  ribs  and  of  skin  stiffeners  and  ribs  were  model 
node  points.  The  node  points  in  turn  defined  the  elements.  The  proce¬ 
dure  for  assigning  area  properties  for  elements,  particularly  for  rods, 
was  detailed  by  both  Heard  (1)  and  Jordan  (33)-  A  brief  summary  is  pre¬ 
sented  in  Appendix  C. 
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5.2  Model  Variations 

The  model  developed  by  Jordan  gave  him  good  response  for  conditions 
where  bending  dominated;  however,  it  provided  no  torsional  stiffness  for 
the  heavy  structural  members.  To  evaluate  damage  and  load  combinations 
producing  significant  torsion,  model  revisions  included  torsional  stiff¬ 
ness  for  the  front  and  rear  spars. 

This  stiffness  was  provided  by  including  rod  elements  along  the 
centerlines  of  the  spars.  These  elements  had  no  axial  load  capacity  but 
did  provide  torsional  resistance.  Multipoint  constraint  equations  deter¬ 
mined  the  rotation  of  each  end  of  a  torsion  rod  by  using  the  lateral  dis¬ 
placements  of  the  nodes  immediately  above  and  below  it.  Figure  7  illus¬ 
trates  that  the  rotation,  6,  of  the  end  of  the  centerline  rod  was 

(5.D 

Although  modeling  philosophies  in  the  previous  efforts  were  essen¬ 
tially  the  same.  Heard  used  membrane  elements  for  the  skin  while  Jordan 
used  shear  panels.  This  study  compared  four  modeling  combinations.  All 
four  models  used  rods  for  skin  stiffeners  and  for  caps  of  spars  ar.d  ribs. 
All  used  shear  panels  for  spar  and  rib  webs.  The  differences  are  pre¬ 
sented  in  Table  I.  Appendix  D  is  a  listing  of  Model  A  and  Appendix  E  is 
a  listing  of  Model  C.  The  additions  for  torsional  resistance  to  convert 
Models  A  and  C  to  Models  8  and  0,  respectively,  are  presented  in  Appendix 
F. 


5-3  Modeling  Initial  Damage 

To  the  extent  possible,  no  special  modeling  techniques  were  applied 


to  initial  damage.  Damaged  shear  panels  or  membranes  were  reduced  in 


Figure  7-  Rod  Elements  for  Spar 
Torsional  Capacity 


TABLE  I 

DESIGNATION  OF  FINITE  ELEMENT  MODELS 


Mode  1 

Designation 

Skin  Elements 

Torsional 

St  i  f fness 

A 

Shear  Panels 

No 

B 

Shear  Panels 

Yes 

C 

Membranes 

No 

D 

Membranes 

Yes 
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thickness  or  were  removed  when  damage  was  severe.  The  rods  representing 
damaged  spars  were  reduced  in  size  to  maintain  an  equivalent  moment  of 
inertia,  as  presented  by  Jordan  (32,  33).  One  exception  to  this  approach 
was  evaluated  for  Test  1.  Damage  to  the  front  spar  extended  halfway  into 
the  web.  Figure  Sa  shows  a  side  view  of  the  damaged  front  spar  and  Fig¬ 
ure  8b  shows  modeling  of  the  undamaged  spar.  The  simpler  modeling  tech¬ 
nique  is  illustrated  in  Figure  8c.  The  web  element  thickness  was  reduced 
to  half  its  undamaged  size.  Rods  representing  spar  caps  were  unmoved  but 
reduced  in  size  to  represent  the  residual  moment  of  inertia.  Figure  8d 
shows  the  more  detailed  approach  used  by  Jordan  (32,  33)  to  model  such 
severe  damage.  The  two  methods  were  compared.  The  specific  changes  made 
to  Models  A  and  C  for  each  test  are  presented  in  Appendix  G.  The  addi¬ 
tional  changes  for  Models  B  and  D  are  included  as  part  of  Appendix  F. 


CHAPTER  VI 


ANALYSIS  AUTOMATION 

6. 1  General 

One  shortcoming  of  the  AFATL  method  cited  in  Chapter  III  was  the 
need  to  examine  voluminous  computer  output.  A  FORTRAN  IV  computer  code, 
entitled  PROSCAN  for  Progressive  Structural  £ol lapse  Analysis,  was  writ¬ 
ten  to  alleviate  the  problem.  PROSCAN  was  written  to  apply  the  method  in 
conjunction  with  NASTRAN  (National  Aeronautics  and  Space  Administration 
Structura 1  Analysis)  to  perform  the  finite  element  analyses.  Figure  9 
illustrates  the  analysis  procedure,  and  Appendix  H  comprises  a  functional 
flow  chart  and  a  listing  of  the  PROSCAN  program. 

In  exchanging  information  between  the  two  computer  programs,  disk 
storage  was  used  exclusively.  All  NASTRAN  output  was  stored  in  punched- 
card  format  in  disk  files.  All  case  control  and  bulk  data  decks  were 
also  stored  on  disks,  and  all  modifications  made  by  PROSCAN  to  the  models 
were  directed  to  those  storage  files. 

6.2  Overstressed  Elements 

The  first  requirement  in  applying  the  AFATL  method  to  finite  element 
results  was  identifying  overstressed  elements.  Because  more  than  one 
limiting  stress  value  was  permissible,  some  common  basis  for  evaluating 
severity  of  stress  had  to  be  established.  The  criterion  selected  was  the 
margin  of  safety  defined  as 
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Figure  9.  Iterative  Analysis  Procedure 
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(6.1) 


Several  elements  within  NASTRAN,  rods  and  shear  panels  included,  re¬ 
turn  an  element  margin  of  safety  as  part  of  the  solution.  For  those  ele¬ 
ments  which  do  not  provide  a  margin  of  safety,  PROSCAN  calculated  one. 
Element  principal  and  maximum  shear  stresses  were  compared  to  analyst- 
provided  limiting  stresses  for  tension,  compression,  and  shear.  A  margin 
of  safety  was  calculated  for  each  type  of  stress  and  the  algebraically 
smallest  value  was  selected  as  the  element  margin  of  safety.  PROSCAN 
then  identified  any  element  with  a  negative  margin  of  safety  as  an  over¬ 
stressed  element. 


6.3  Failed  Elements 

PROSCAN  applied  the  next  step  of  the  AFATL  method,  grouping  of  over¬ 
stressed  elements,  by  node  matching.  The  node  numbers  of  each  overstress¬ 
ed  element  were  compared  to  those  of  every  other  overstressed  element. 
PROSCAN  designated  any  continuous  linkage  of  those  elements  as  a  group, 
then  selected  the  most  severely  stressed  element  from  the  group. 

The  element  margin  of  safety  again  was  the  basis  for  decisions.  PRO¬ 
SCAN  selected  the  element  of  the  group  with  the  most  negative  margin  of 
safety.  That  element  became  a  failed  element.  The  process  of  grouping 
and  failing  elements  continued  until  all  overstressed  elements  were  con- 
s i dered. 

PROSCAN  did  not  actually  remove  a  failed  element  from  the  model. 
Instead,  PROSCAN  assigned  property  values  to  the  element  which  effective¬ 
ly  eliminated  its  contribution  to  the  structure.  The  failed  areas  and 
moduli  of  elasticity  and  shear  were  orders  of  magnitude  below  nominal 
values  for  unfailed  elements. 
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6.4  Propagation  of  Damage 

The  failing  of  an  element  represented  propagation  of  the  damage,  and 
as  a  consequence,  the  borders  of  the  damage  expanded.  Additional  ele¬ 
ments  had  to  be  identified  as  bordering  the  new  damage  so  they  could  be 
assigned  reduced  limiting  stresses.  Again  a  node  matching  scheme  was  em¬ 
ployed.  Each  element  which  had  at  least  one  node  in  common  wi th  a  newly 
failed  element  was  examined.  If  it  had  not  already  failed  itself  or  had 
not  already  bordered  damage,  lower  limiting  stresses  replaced  those  pre¬ 
viously  used.  The  lowering  of  limiting  stresses  accounted  for  the  possi¬ 
ble  presence  of  crack  tip  stress  concentrations  as  introduced  in  section 
2.5. 


6.5  Adjustment  of  Load 

PROSCAN  had  the  capability  of  applying  a  new  load  to  the  model  with 
each  iteration.  That  capability  was  used  in  this  study  as  explained  be¬ 
low. 

If  no  element  failed  on  a  particular  iteration,  the  load  was  increas¬ 
ed  for  the  next  NASTRAN  analysis.  This  would  occur  until  the  structure 
sustained  some  maximum  user-specified  load  without  further  element  fail¬ 
ure.  Conversely,  if  an  element  failed  on  a  particular  iteration,  PROSCAN 
reduced  the  load  for  the  next  NASTRAN  analysis.  The  purpose  was  to  deter¬ 
mine  the  structure's  ability  to  carry  a  lesser  load  after  further  weaken¬ 
ing  by  the  failed  element.  Reducing  the  load  every  time  an  element  fail¬ 
ed  continued  until  the  structure  could  not  sustain  a  minimum  load  without 
further  fa i lure. 

The  analyst  provided  a  sequence  of  loads  to  be  applied,  from  minimum 
to  maximum,  as  part  of  the  PROSCAN  input  data.  PROSCAN  then  made  the 
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appropriate  changes  to  the  NASTRAN  case  control  deck  to  reflect  the 
structure's  performance  on  the  previous  iteration.  In  addition  to  chang¬ 
ing  the  load  identification  number,  PROSCAN  could  assign  new  single  point 
and  multipoint  constraint  sets  and  identify  new  labels  to  correspond  to 
each  new  load. 

PROSCAN  automated  the  entire  application  of  the  AFATL  method.  This 
began  with  initial  viewing  of  NASTRAN  output  and  finished  by  establishing 
new  files  containing  modified  case  control  and  bulk  data  decks. 
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CHAPTER  VI  I 


DETERMINATION  OF  LIMITING  STRESSES 
7.1  Need  for  Limiting  Stresses 

Repeated  reference  has  been  made  to  limiting  stresses.  It  is  appro¬ 
priate  to  address  in  more  detail  the  specifics  of  allowable  stress  levels. 
Heard  (1)  used  two  limiting  stress  criteria:  ultimate  strength  for  ele¬ 
ments  away  from  damage,  and  yield  strength  for  elements  bordering  damage. 
This  study  attempted  to  define  more  precisely  the  levels  of  stress  which 
should  cause  failure  in  the  model. 

Ultimate  strength  remained  the  basic  criterion  for  defining  failure; 
but  most  elements,  even  those  away  from  damage,  were  assigned  limiting 
stresses  lower  than  ultimate  strength.  Consider  that  a  relatively  large 
element  returned  a  computed  stress  representative  of  a  large  structural 
region.  This  representative  stress  was  unavoidably  lower  than  the  high 
stress  within  the  region  which  would  cause  failure  in  the  actual  member. 

It  was  necessary  then  to  estimate  the  effects  of  the  representative 
st  resses  by  using  some  value  of  1  imi  t  ing  st  ress  1  ower  than  ultimate  strength . 

Appropriate  limiting  stresses  also  estimated  the  nonlinear  behavior 
experienced  through  the  buckling  of  skin  panels.  Since  the  finite  element 
analyses  assumed  linearly  elastic  behavior,  the  ability  to  compensate  for 
skin  panel  buckling  was  incorporated  to  enhance  results.  PROSCAN  had  the 
ability  to  incorporate  both  the  low  stresses  causing  buckling  and  the  re¬ 
duced  stiffnesses  subsequent  to  buckling. 
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7.2  Limiting  Stresses  for  Rod  Elements 

Spar  and  rib  sections  were  each  represented  by  three  elements.  A 
shear  panel  represented  the  web.  One  rod  element  represented  the  upper 
spar  cap  and  another  rod  represented  the  lower  spar  cap.  The  rods  were 
sized  and  spaced  to  maintain  the  moment  of  inertia  about  the  section's 
neutral  axis  and  to  return  outer  fiber  stresses. 

The  stress  value  obtained  for  rod  elements  was  the  average  of  the 
stresses  at  each  end  of  the  rod.  Because  rod  elements  in  the  spars  were 
relatively  long,  the  average  stress  could  be  substantially  less  than  the 
maximum  stress.  A  procedure  to  obtain  limiting  stresses  for  a  similarly 
modeled  doubly  symmetric  cantilevered  beam  served  as  a  foundation  for 
developing  limiting  stresses  for  the  wing  model.  Figure  10  shows  such  a 
beam  with  top  rod  elements  numbered  and  top  nodes  lettered. 


Lo*d  P 


Figure  10.  Cantilevered  Beam  Model 


Using  rod  3  for  illustration,  the  maximum  stress  from  the  applied 


load  occurred  at  node  E,  but  the  stress  obtained  was  the  average  of 


stresses  at  nodes  0  and  E.  Designating  U  as  the  length  of  any  member  j, 
the  average  stress  in  rod  3  was 


P(L,  ♦  L2  ♦  1  !.,)» 


avg  I 

and  the  stress  at  node  E  was 


(7.1) 


P(L,  +  L2  +  L3)y 


max 


(7.2) 


where  I  was  the  section  moment  of  inertia.  Designating  a  limiting  stress 


factor,  F. ,  as  the  ratio  of  a  to  a  , 
i  avg  max 

F  .  Li  *  4  *  1 4 

3  -  L,  *  L2  *  4 


(7.3) 


In  general  terms, 


1  M 

i  4 

j=i  J 


(7.4) 


for  the  single  point  load  shown.  The  appropriate  limiting  stress,  , 

i 

was 


°l  “  F-  CT  H 

L.  i  ul t . 


(7.5) 


where  a  .  was  the  ultimate  strength  for  the  member  i. 

To  extend  this  approach  for  calculating  stress  factors  to  the  finite 
element  model  of  the  wing,  the  wing  itself  was  idealized  as  a  straight 
cantilevered  beam.  The  front  spar  dimensions  were  used  for  section 
lengths  as  shown  in  Figure  11.  Upper  surface  element  numbers  are  below 
each  rod  and  corresponding  node  numbers  are  above  each  node. 
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Figure  11.  Cantilevered  Front  Spar  Idealization 


Factors  to  reduce  material  ultimate  strength  for  elements  inboard  of 
the  load  were  calculated  as  shown  in  the  previous  example.  Elements  be¬ 
tween  the  load  and  the  wing  tip  used  the  same  factor  as  the  elements  im¬ 
mediately  inboard  of  the  load.  Table  II  shows  the  limiting  stress  fac¬ 
tors  for  the  front  spar  rods  on  the  upper  wing  surface. 


TABLE  11 

LIMITING  STRESS  FACTORS  FOR  ROD  ELEMENTS 


Rod  No.:  641  613  575  529  475  431  393  355 

F.:  0.350  0.937  0.921  0.855  0.738  0.735  0.500  0.500 


Each  rod  representing  a  skin  stiffener  was  approximately  parallel  to 
the  spars  and  was  assigned  the  same  factor  as  its  corresponding  spar  ele¬ 
ment.  Each  rib  was  approximately  perpendicular  to  the  spars.  Each  rod 
in  a  rib  was  assigned  the  factor  of  the  spar  rod  immediately  inboard  of 
the  spar-rib  intersection. 


A  single  concentrated  load  was  used  for  analysis  to  correspond  to 
the  actual  loading  applied  in  the  laboratory  test  program.  However,  an 
aerodynamic  load  could  be  represented  by  any  approximation  acceptable  to 
the  analyst.  Although  the  mathematical  expression  for  F.  would  be  more 
complex,  the  same  approach  to  factoring  for  limiting  stresses  in  rod  ele¬ 
ments  could  be  applied. 


7.3  Limiting  Stresses  for  Web  Elements 


Shear  panel  elements  represented  the  webs  of  spars  and  ribs.  The 
limiting  stress  for  shear  was  determined  by  comparing  the  average  shear¬ 
ing  stress  in  the  web  to  the  maximum  shearing  stress  in  the  web.  If  V 
were  designated  as  the  shearing  force  in  the  cantilevered  beam  discussed 
in  the  previous  section,  the  web  element  yielded  a  shearing  stress  of 


T 

avg 


V 

2yt 


(7.6) 


where  t  was  the  web  thickness.  The  maximum  shearing  stress  in  the  sec¬ 
tion  was 


T 

max 


(7.7) 


The  limiting  stress  factor,  F,  was  the  value  of  x  divided  by  x  ,  so 
3  avg  '  max 

the  limiting  shearing  stress,  x^,  was 


ult 


(7.8) 


Calculations  for  typical  spar  cross  sections  showed  F  =  0.35  to  be  a 
representative  value.  This  value  was  applied  to  all  spar  and  rib  web 


elements. 
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Proportions  of  spar  and  rib  sections  indicated  web  crippling  was  un¬ 
likely;  therefore,  no  reductions  in  limiting  stresses  were  developed  for 
buckling  of  undamaged  web  members.  If  initial  damage  to  the  structure 
introduced  a  potential  for  buckling,  residual  member  proportions  dictated 
the  appropriate  reductions. 

7.4  Limiting  Stresses  for  Skin  Elements 

Skin  panels,  unlike  spar  and  rib  webs,  were  susceptible  to  buckling. 
Additionally,  the  skin  could  tear  along  rivet  lines  or  rivets  themselves 
could  fail.  Whether  a  panel  buckled  in  shear  or  in  compression,  or  fail¬ 
ed  along  a  rivet  line,  the  result  was  a  reduction  in  stiffness  of  the 
panel.  Because  of  the  similar  change  in  behavior,  panels  were  divided 
into  either  ore-buckling  or  post-buckling  categories  even  though  rivet 
line  failure  was 'not  a  buckling  phenomenon. 

Each  skin  panel  on  the  wing  had  slightly  different  geometric  proper¬ 
ties  which  gave  each  slightly  different  pre-  and  post-buckling  character¬ 
istics.  Panel  106,  forward  of  the  front  spar  on  the  lower  wing  surface, 
was  typical  of  most  skin  panels  and  was  used  to  determine  approximate 
values  for  all  panels.  Figure  12  shows  its  location  in  the  model.  The 
assembly  used  for  calculations  included  panel  106,  a  stiffener  attached 
to  each  long  side,  and  a  rib  attached  to  each  short  side.  Averaging  the 
lengths  of  the  two  long  sides  and  the  two  short  sides  gave  a  rectangular 
shape  for  calculations. 

To  determine  pre-buckling  limits,  compression  perpendicular  to  the 
long  sides,  compression  perpendicular  to  the  short  sides,  and  a  corner 
force  producing  shear  were  all  evaluated  separately.  Calculations  were 
determined  according  to  Peery  (34,  Chapters  14  and  15).  The  average 
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stress  causing  buckling  in  each  case  was  divided  by  the  material  ultimate 
strength  to  determine  the  limiting  stress  factors,  F.  For  the  two  com¬ 
pression  conditions,  the  more  conservative  value  was  used. 


Figure  12.  Typical  Panel  for  Buckling  Limits 


Limits  for  tension  were  obtained  by  calculating  rivet  and  skin 
strengths  along  a  conservative  rivet  line.  Limiting  loads,  determined 
according  to  Peery  (3h,  Chapter  12)  and  Bruhn  (35,  Chapter  01),  were 
divided  by  the  ultimate  load,  the  load  causing  an  average  stress  in  the 
panel  equal  to  the  ultimate  strength.  The  result  was  the  limiting  stress 
factor.  Skin  failure  was  compared  to  rivet  failure,  and  the  more  conserv 
ative  value  was  used. 

For  oost-buckl ing  behavior,  the  limiting  stresses  were  returned  to 
ultimate  strength,  but  the  elastic  and  shear  moduli  were  reduced  to 
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account  for  the  reduction  in  stiffness  following  buckling.  For  compres¬ 
sion  buckling,  bilinear  behavior  was  assumed.  The  load  versus  displacement 
curve  of  Figure  13a  assumed  linear  behavior  prior  to  buckling,  then  linear 
behavior  from  buckling  to  an  ultimate  strength  failure.  The  slope  of  the 
pre-buckling  portion  of  the  curve,  S  ,  was  divided  into  a  secant  slope, 
from  the  origin  to  failure.  The  result  was  a  reduction  factor,  M,  for  the 
elastic  modulus.  The  same  process  applied  to  Figure  1 3b  produced  a  reduc¬ 
tion  factor  for  the  modulus  of  shear.  Table  III  summarizes  the  results 
for  ski n  buckl i ng. 


TABLE  I  I  I 

LIMITING  FACTORS  FOR  SKIN  ELEMENTS 


Behavior 

Stress  Factor 

F 

Modulus  Factor 
M 

Compress i on , 
Pre-buckl i ng 

0.17 

1 .00 

Shear , 

Pre-buck  1 i ng 

0.  32 

1  .00 

Tens i on , 

Pre-buckl i ng 

0.55 

o 

o 

Compression, 
Post-buck! i ng 

1.00 

0.76 

Shear , 

Post-buckl i ng 

1.00 

0.53 

Tens i on , 

Post-buckl i ng 

1.00 

0.76 

Ditplacement 


38 


Although  the  assumption  of  linear  behavior  from  buckling  to  ultimate 
failure  was  not  correct,  it  was  a  conservative  representation  of  the 
rather  brittle  material  behavior  observed  in  the  laboratory.  The  result 
economically  approximated  the  loss  in  stiffness  suffered  by  the  structure 
from  skin  buckling  and  rivet  line  failure. 

7.5  Damage  Propagation 

No  attempt  was  made  to  model  ragged  edges  around  initial  damage  nor 
to  reduce  element  size  in  areas  of  propagating  cracks.  The  large  ele¬ 
ments  then  tended  to  mask  the  stress  concentrations  around  cracks  and  pro¬ 
duced  a  model  significantly  more  resistant  to  progressive  collapse  than 
the  structure  being  represented. 

Conventionally,  the  nominal  stress  in  a  cracked  member  would  have 
been  multiplied  by  a  stress  concentration  factor,  K.  Its  value  would  have 
been  larger  than  1.0  and  based  upon  crack  length  and  crack  tip  severity. 
The  increased  value  for  stress  at  the  crack  tip  would  then  have  been  com¬ 
pared  to  an  allowable  stress  for  the  member.  PROSCAN  used  an  inverse  ap¬ 
proach.  Rather  than  increase  the  nominal  stress  returned  by  an  element, 
the  allowable  stress  was  decreased  by  a  factor  F,  where  F  essentially  was 
the  inverse  of  K.  This  further  reduction  of  limiting  stresses  compensat¬ 
ed  for  the  absence  of  increased  modeling  detail  around  damage. 

Any  cracks  occurring  were  assumed  to  originate  at  the  initial  damage 
or  in  subsequently  failed  elements.  The  further  reductions  in  limiting 
stresses  applied  therefore  only  to  unfailed  elements  bordering  either 
initial  damage  or  failed  elements. 

Separate  reduction  factors  were  determined  for  tension  and  for  shear. 
Because  cracks  were  assumed  not  to  propagate  in  compression,  no  further 
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reduction  applied  to  limiting  compressive  stresses.  This  portion  of  the 
investigation  was  patterned  after  similar  crack  propagation  studies  by 
Sih  and  Hartranft  (28). 

Two  square  plate  models,  one  loaded  in  tension  and  the  other  in 
shear,  provided  information  for  reduction  factors.  Model  detail  ranged 
from  two  elements  along  a  side  to  thirty-two  elements  along  a  side. 

A  crack  initiated  at  the  center  of  one  edge  propagated  through  the 
plate  during  sequential  analyses.  Loads  remained  constant  through  all 
iterations.  Crack  propagation  was  represented  by  creating  a  new  node  be¬ 
side  the  tip  of  the  crack,  thus  extending  the  crack  to  the  next  node. 
Figure  14  illustrates  the  procedure,  exaggerated  in  scale,  on  a  model 
using  four  elements  per  side.  Figure  14b  shows  node  m  at  the  tip  of  the 
crack.  The  creation  of  node  z  extended  the  crack  tip  to  node  n  in  Fiqure 
14c. 

Figure  1 4b  shows  the  plate  cracked  one-quarter  of  the  way  through 
its  width.  Stresses  in  the  four  elements  connected  to  the  node  at  the 
crack  tip,  those  indicated  by  X's,  were  averaged  and  then  di vided  into  the 
average  stress  in  the  uncracked  plate.  The  result  was  the  limiting  stress 
factor  for  the  plate  cracked  through  one-quarter  of  its  width.  The  same 
procedure  applied  to  the  plate  in  Figure  14c  produced  the  factor  for  the 
plate  cracked  halfway  through  its  width.  Figure  15  shows  variation  of 
the  factor  as  a  function  of  model  detail  and  crack  length.  F  represents 
tension  loading  and  F^  represents  shear  loading. 

All  curves  in  Figure  15  appeared  to  approach  zero  slope  as  element 
size  reduced.  The  values  of  F^  and  F^  selected  for  this  study  were  for 
a  plate  cracked  one-eighth  of  the  way  through  its  width  in  a  model  with 
32  elements  along  a  side.  For  most  components  of  the  F-84F  wing,  this 


Figure  14.  Crack  Investigation  Models 


Number  of  Elements  per  Side 
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represented  a  crack  less  than  one  inch  in  length  in  the  longer  side  of 
the  component.  The  corresponding  stress  concentration  factor  was  K  =  2.1. 
For  skin  panel  106,  it  represented  a  crack  0.725  inches  long  with  a  crack 
tip  radius  of  0.30  inches  calculated  according  to  Seely  and  Smith  (36, 
Chapter  12).  Thus  the  assumed  crack  around  damage  was  relatively  mild 
and  was  therefore  conservative. 

New  limitinq  stresses  for  an  element  bordering  damage  were  the  pro¬ 
duct  of  the  appropriate  factor,  F^  or  F  or  F  =  1.0  for  compression, 
and  the  element's  previous  limiting  stresses.  Limiting  stress  factors 
were,  in  this  manner,  cumulative.  The  exception  was  unbuckled  skin  pan¬ 
els.  Their  limiting  stresses  were  not  reduced  to  reflect  cracks  until 
after  buckling  stresses  were  exceeded. 


CHAPTER  VIII 


COMPARISON  OF  RESULTS 

8.1  General 

Damage  and  load  conditions  for  Tests  2C  and  3B  were  analyzed  using 
all  four  models  for  each  test.  Measured  rotations  of  wing  spar  roots 
from  laboratory  data  were  enforced  in  the  analyses.  Examination  of  those 
analyses  showed  the  addition  of  torsional  rod  elements  to  the  spars  made 
little  difference  in  results.  The  performance  of  Model  B  was  very  simi¬ 
lar  to  that  of  Model  A,  and  the  results  from  Model  D  were  almost  identi¬ 
cal  to  those  from  Model  C.  Apparent  reasons  for  the  similarities  are 
presented  in  the  next  section. 

Further  examination  of  the  analytical  results  revealed  unexpected 
stress  distributions  in  and  near  the  wing  spar  roots.  The  enforced  rota¬ 
tions  of  wing  spar  roots,  although  developed  from  experimental  measure¬ 
ments,  did  not  produce  purely  rigid  body  motions  for  reasons  explained  in 
section  8.6.  Consequently,  the  original  analytical  representations  did 
not  match  closely  enough  the  laboratory  conditions  of  the  experimental 
test  program. 

A  second  set  of  analyses  was  performed  using  zero  support  rotations 
for  stress  determination  and  enforced  rotations  for  checking  displace¬ 
ments.  Tests  2C  and  were  analyzed  using  Models  A  and  C  for  each  test. 
Model  C  described  the  collapse  phenomenon  more  closely  than  Model  A  as 
explained  in  section  8.1*.  Therefore,  Model  C  was  next  compared  to  Model 
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0,  the  same  model  with  the  addition  of  torsional  rod  elements.  The  lack 
of  significant  different  between  Models  C  and  D  confirmed  the  minimal  in¬ 
fluence  of  the  torsional  rod  elements.  Model  B,  therefore,  was  not  ana¬ 
lyzed  further  because  it  would  produce  essentially  the  same  results  as 
Model  A.  Even  though  torsional  rod  elements  were  not  significantly  af¬ 
fecting  results,  Model  D  was  selected  for  the  comparison  of  initial  dam¬ 
age  modeling  since  its  torsional  capability  could  provide  greater  lati¬ 
tude  for  an  analyst  to  adjust  model  stiffness. 

Model  D  was  used  to  evaluate  the  two  approaches  to  modeling  damage 
for  Test  1  described  in  section  5-3-  The  simpler  method  of  modeling  por¬ 
trayed  more  accurately  the  pattern  of  failure  as  explained  in  section 
8.5.  The  simpler  method  of  modeling  the  damage  was  then  applied  to  Model 
A  for  a  final  analysis  of  Test  1. 

8.2  Comparison  of  Failure  Loads 

A  close  correlation  of  analytically  predicted  failure  loads  with  ex¬ 
perimentally  measured  failure  loads  would  be  a  desirable  result  of  evalu¬ 
ating  the  AFATL  method.  Table  IV  summarizes  the  failure  load  results. 

The  models  ranged  from  5  percent  to  85  percent  stronger  than  the  actual 
structure.  Note  that  for  Test  2  no  experimental  failure  load  was  deter¬ 
mined;  therefore,  conclusions  about  Test  2  are  judgmental. 

Mode!  A  gave  the  closest  approximation  for  Test  3  and  may  have  given 
a  close  approximation  for  Test  2.  However,  for  reasons  discussed  in  sec¬ 
tion  8.4,  Model  A  was  not  considered  the  best  model.  Model  D  was  more 
conservative  than  Model  A  in  estimating  wing  strength.  Although  Model 
0's  predicted  strength  for  Test  2  was  clearly  less  conservative  than  for 
Test  3i  the  results  may  have  been  acceptably  consistent.  Both  approaches 
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for  modeling  Test  1  damage  gave  excessive  predicted  strengths;  however, 
section  8.5  discusses  how  those  figures  might  be  improved. 


TABLE  IV 

SUMMARY  OF  FAILURE  LOADS 


Test 

Fai 1 ure  Loads  (k i ps) 

Laboratory 

Mode  1  A 

Model  C 

Mode  1  D 

1 

12.0 

22 

— 

22 

2 

15-0* 

18 

20 

20 

3 

12.4 

13 

19 

19 

Largest  load  applied;  no  failure  load  determined. 


Models  C  and  0  showed  no  differences  in  failure  loads  and  very  lit¬ 
tle  difference  in  the  sequences  of  element  failure.  There  are  two  appar¬ 
ent  reasons  for  the  similarity.  The  first  is  that  the  torsional  capaci¬ 
ties  of  the  spars  were  probably  underestimated  when  the  torsional  rod 
elements  were  sized.  Second,  bending  was  the  dominant  behavior  of  the 
F-84F  wing  even  under  extreme  conditions  such  as  those  of  Test  2. 

8.3  Load-Iteration  History 

The  AFATL  method,  as  applied  by  PROSCAN,  caused  loads  to  vary  from 
iteration  to  iteration.  Figure  16  depicts  the  variation  of  load  with  re¬ 
spect  to  iteration  for  the  first  35  cycles  for  Models  A  and  D.  The  analy¬ 
tical  data  for  any  given  load  level  were  taken  from  the  last  cycle  in  which 
that  load  was  applied  before  the  model  experienced  a  higher  load.  For 
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Figure  16.  (Continued) 


example,  the  analytical  data  for  the  Model  A  analysis  of  Test  2C  at  16 
kips  applied  load  came  from  iteration  No.  8.  As  shown  in  Figure  16c, 
that  was  not  the  first  application  of  a  16-kip  load,  but  it  was  the  last 
iteration  before  a  higher  load,  18  kips,  was  applied. 

That  procedure  for  selecting  which  iterations  to  use  for  data  com¬ 
parison  occasionally  led  to  gaps  of  several  iterations  between  successive 
data-produc i ng  loads.  Again  as  an  example.  Figure  1 6 f  shows  13  iterations 
elapsed  between  the  1 4-kip  and  16-kip  loads  for  the  Model  D  analysis  of 
Test  38.  During  those  cycles,  six  elements  failed.  This  characteristic 
of  the  procedure  accounted  for  the  occasional  sharp  discontinuities  in 
the  plots  of  data. 

Figure  16  also  emphasizes  the  need  for  caution  in  setting  the  mini¬ 
mum  load  to  be  investigated.  PROSCAN  permitted  the  load  to  drop  consid¬ 
erably  during  a  series  of  element  failures,  then  again  rise  to  a  high 
level.  Figure  16a  shows  how  the  load  dropped  from  22  kips  down  to  12 
kips  before  again  climbing  back  up  to  22  kips.  Making  the  minimum  allow¬ 
able  load  too  large  could  result  in  a  premature  indication  of  structural 
failure.  It  could  occur  during  such  a  series  of  element  failures  when, 
in  fact,  the  structure  still  possessed  the  capacity  for  loads  well  above 
the  mini  mum  I  eve  1 . 

8.4  Internal  Load  Paths 

The  most  demanding  test  of  the  models  was  how  realistically  they 
transferred  the  loads  internally  through  the  wing  structure  and  into  the 
supports.  Figure  26  (Appendix  l)  compares  the  vertical  support  reactions 
for  experimental  and  analytical  results.  Figures  27  through  29  (also 
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Appendix  I)  compare  variations  of  strain  at  representative  points  on  the 
wing  with  respect  to  applied  load. 

Examination  of  Figure  26  through  29  showed  that  neither  Model  A  nor 
Model  D  transferred  the  load  from  the  loaded  spar  to  the  unloaded  spar 
as  quickly  as  the  actual  wing  did.  Additionally,  neither  model  trans¬ 
ferred  as  much  of  the  load  from  spar  to  spar  as  the  wing  did. 

The  most  important  indication  for  this  study  of  how  realistically 
the  models  transferred  the  loads  internally  came  from  Figures  17  through 
19-  They  depict  the  buckled  and  failed  elements  in  Models  A,  C,  and  D 
at  their  respective  failure  loads.  For  Test  3B,  Figure  19,  Model  A  did 
not  indicate  the  nature  of  the  failure  as  observed  in  the  experimental 
test  program;  however,  Models  C  and  D  did  match  closely  the  laboratory 
observations.  For  Test  2C,  Figure  18,  no  failure  occurred  in  the  experi¬ 
mental  program,  but  Models  C  and  D  predicted  a  plausible  failure.  Model 
A,  however,  predicted  failure  of  the  front  spar  at  one  of  its  strongest 
sections.  For  Test  1,  Figure  17.  Model  D  matched  the  laboratory  failure 
pattern  very  closely  using  the  simple  modeling  of  initial  damage.  Model 
A,  however,  indicated  failure  of  the  undamaged  rear  spar.  All  models 
indicated  more  overstressing  of  skin  elements  near  the  wing  spar  roots 
than  was  observed  on  the  actual  wing.  A  complete  summary  of  results  for 
the  first  35  iterations  of  the  principal  series  is  presented  in  Appendix 
J. 

8.5  Comparison  of  Damage  Modeling 

Section  5-3  introduced  two  approaches  for  modeling  Test  1  damage. 
Both  approaches  predicted  the  same  failure  load,  but  Figure  17  illus- 
strates  that  there  were  significant  differences  in  which  elements  failed. 
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(b)  Model  C  (without  torsional  stillness  rods) 

Figure  18.  (Continued) 
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failed  Membrane  Elements 
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(b)  Model  C  (wit  ho  til  iofsion.il  stillness  rods) 

Figure  19.  (Continued) 
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Of  the  two  approaches  using  Model  D,  the  simpler  approach  depicted  more 
accurately  the  failure  as  observed  in  the  laboratory.  The  simple  ap¬ 
proach  applied  to  Model  A  did  not  produce  an  accurate  failure  pattern; 
however,  all  tests  indicated  that  Model  A  was  less  suitable  for  predict¬ 
ing  the  pattern  of  failure. 

Although  the  Model  0  results  suggested  a  preference  for  the  simpler 
technique,  caution  is  advised  before  reaching  a  firm  conclusion.  In  all 
models  for  all  tests,  the  webs  of  spars  and  ribs  were  represented  by 
shear  panels.  Consequently,  the  front  spar  in  Test  1  could  not  fail  at 
the  damage  until  shear  limits  were  exceeded.  The  experimental  program 
showed  the  damaged  front  spar  web  in  Test  1  failed  in  bending  tension,  a 
failure  mode  the  shear  panel  could  not  predict.  For  cases  of  initial 
damage  where  all  or  most  of  a  spar  cap  or  rib  cap  would  be  removed,  the 
web  should  probably  be  modeled  by  a  membrane  element.  Although  the  mem¬ 
brane  element  would  be  stiffer  than  the  shear  panel,  it  would  be  direct¬ 
ly  sensitive  to  limiting  tensile  and  compressive  stresses  as  well  as  to 
shear  limits.  Such  a  recommendation  applied  to  this  study  may  have  ap¬ 
preciably  reduced  the  predicted  failure  load  for  Test  1,  and  it  may  have 
altered  the  apparent  value  of  simple  modeling  over  the  more  detailed 
representat ion. 

8.6  Rotation  of  Wing  Spar  Roots 

Specific  values  of  displacement  were  of  interest  in  this  study  as  an 
additional  means  of  comparing  analytical  results  to  laboratory  data.  To 
obtain  more  accurate  displacement  values  from  the  analytical  method,  rota¬ 
tions  of  wing  spar  roots  were  measured  in  the  experimental  test  program 
and  enforced  in  the  analytical  models.  However,  for  most  applications  of 
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the  AFATL  method,  rotations  at  structure  supports  would  not  be  known. 
Additionally,  the  precise  displacements  of  the  structure  probably  would 
be  unimportant.  The  displaced  shape  of  the  structure,  which  might  be 
used  to  modify  loading  for  each  iteration,  was  available  from  the  analy¬ 
ses  using  zero  support  rotations. 

The  enforced  rotations  were  derived  from  experimental  data.  In 
translating  laboratory  measurements  into  single  point  constraints  for 
NASTRAN,  an  assumption  was  made.  It  was  assumed  that  the  center  of  rota¬ 
tion  for  each  spar  was  the  point  midway  between  the  two  pins  securing  the 
spar  in  its  support  structure.  In  fact,  any  point  between  those  two  pins 
could  have  been  the  center  of  rotation,  and  the  center  could  have  changed 
as  loading  progressed.  The  assumption  almost  certainly  contributed 
to  the  introduction  of  erroneous  stresses  into  the  models  during  the 
first  set  of  analyses. 

Ano'her  likely  contributor  to  those  stresses  was  the  manner  in  which 
some  of  the  multipoint  constraint  equations  for  the  models  were  written. 

A  spar  root  was  modeled  by  a  shear  panel  and  two  rod  elements,  a  configur¬ 
ation  that  gave  the  desired  resistance  to  bending  but  provided  no  lateral 
restraint.  The  necessary  lateral  restraint  was  provided  by  multipoint 
constraints  to  keep  each  root  section  in  line  with  its  adjacent  spar  sec¬ 
tion. 

The  multipoint  constraint  equations  for  the  original  model  were  for¬ 
mulated  not  in  a  general  manner,  but  with  an  implied  assumpt ion  that  there 
was  no  displacement  of  the  wing  spar  root  nodes.  Thus  any  attempt  to  en¬ 
force  the  measured  rotations  violated  that  assumption.  The  result  was 


erroneous  stresses  near  the  base  of  the  wing. 
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Of  the  two  sources  of  error  identified,  the  multipoint  constraint 
equations  could  be  easily  corrected.  The  problem  of  precisely  measuring 
wing  spar  root  rotations  cannot  be  solved  without  sophisticated  measuring 
equipment.  The  benefits  gained  from  precise  measurements,  however,  would 
not  begin  to  justify  the  added  expense  for  normal  applications  of  the 
method. 


8.7  Deflections 

Deflections  were  measured  in  the  experimental  test  program  and  were 
compared  to  analytical  results.  Figures  30  through  32  (Appendix  I)  pre¬ 
sent  single-point  deflection  data  for  Models  A,  C,  and  D,  and  for  the 
test  program.  Although  Model  D  was  selected  as  the  best  model  because 
of  its  ability  to  predict  the  failure  most  realistically,  Model  A  was 
superior  for  predicting  displacement  values.  For  general  deflected 
shape,  however,  there  was  little  difference  between  Models  A  and  D. 

Figure  20  compares  Test  3B  profiles  of  the  front  and  rear  spars  for 
Models  A  and  D,  and  for  the  actual  wing  at  their  respective  failure 
loads.  Both  models  presented  essentially  the  same  deflected  shape  which 
differed  only  slightly  from  measured  results. 

As  mentioned  in  the  previous  section,  deflections  are  not  envisioned 
as  a  critical  factor  in  the  routine  application  of  the  AFATL  method.  Even 
if  deflected  shape  were  important,  Models  A  and  D  both  returned  approxi¬ 
mately  the  same  results.  If  specific  values  of  displacement  were  to  be¬ 
come  the  overriding  concern  in  a  specialized  application,  Model  A  would 
appear  to  be  the  better  model.  Otherwise,  Model  D  provided  reasonable 
accuracy  for  deflected  shape. 
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Figure  20.  Failure  Load  Profiles  for  Test  3B 


CHAPTER  IX 


SUMMARY  AND  CONCLUSIONS 

The  principal  goal  of  this  study  was  to  evaluate  the  suitability  of 
the  AFATL  method  for  predicting  progressive  collapse  in  complex  struc¬ 
tures.  Suitability  was  to  be  investigated  by  determining  supportable 
limiting  stress  values,  selecting  a  good  combination  of  finite  elements 
for  modeling,  and  comparing  analytical  to  experimental  results. 

Limiting  stresses  used  in  this  study  were  a  direct  application  of 
classical  theory.  Consequently,  any  skilled  analyst  could  apply  the  con¬ 
cepts  to  any  structure.  A  conscious  effort  was  made  to  found  the  work  in 
commonly  known  principals  of  materials  behavior  and  to  avoid  the  structure- 
dependency  associated  with  empirical  formulations. 

All  models  evaluated  used  axial  rod  and  shear  panel  element  combina¬ 
tions  to  represent  spars  and  ribs.  Model  D  used  membrane  elements  to 
model  aircraft  skin  and  rod  elements  to  model  skin  stiffeners.  Addition¬ 
ally,  it  had  torsional  rod  elements  along  the  spar  centerlines.  Model  A 
used  shear  panels  and  thickened  rod  elements  to  represent  aircraft  skin 
and  skin  stiffeners.  Model  A  had  no  torsional  rod  elements. 

Models  A  and  D  both  overestimated  the  residual  strength  of  the  dam¬ 
aged  structures.  For  the  purposes  of  this  study,  those  results  were  con¬ 
servative.  A  deficiency  observed  in  the  study  was  the  lack  of  consistency 
in  the  degree  to  which  residual  strength  was  overestimated.  However,  all 
estimates  were  within  a  factor  of  two  of  the  experimental  results. 
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Model  A  provided  better  deflection  estimates  than  Model  D.  Using 
shear  panel  elements  for  the  aircraft  skin  made  Model  A  more  difficult  to 
prepare  than  the  models  using  membrane  elements.  As  explained  in  Appen¬ 
dix  C,  the  use  of  shear  panel  elements  required  additional  calculations 
for  modifying  rod  element  sizes  to  represent  the  membrane  capacity  of  the 
skin.  However,  once  Model  A  was  developed,  it  was  less  expensive  to  use 
than  models  with  membrane  skin  elements.  For  applications  where  struc¬ 
ture  displacements  are  of  primary  concern,  Model  A  would  provide  better 
resul ts. 

Model  0  described  more  accurately  the  actual  pattern  of  failure  lead¬ 
ing  to  structural  collapse.  Using  membrane  elements  for  the  aircraft  skin 
made  Model  D  a  simpler  model  to  prepare  as  described  in  Appendix  C.  The 
membrane  elements  also  gave  a  better  qualitative  representation  cf  skin 
panel  behavior.  For  applications  where  the  failure  pattern  is  of  princi¬ 
pal  interest,  Model  D  would  provide  better  results. 

PROSCAN  was  developed  as  a  convenience  to  automate  the  application 
of  the  AFATL  method.  It  proved  to  be  more  of  a  necessity  than  a  conven¬ 
ience  in  processing  the  volumes  of  data  generated  by  many  iterative  finite 
element  analyses.  Additionally,  it  provided  flexibility  in  the  selection 
of  loading  sequences  and  in  the  application  of  limiting  stresses  for  ele¬ 
ments  . 

The  combination  of  automation,  modeling  techniques,  and  limiting 
stresses  applied  to  the  AFATL  method  produced  a  useful  estimating  tool 
for  predicting  progressive  collapse  in  complex  structures  such  as  the 
F-84F  aircraft  wing.  The  F-8AF  wing  is  a  semi -monocoque  structure  with  a 
heavy  two-spar  skeletal  frame.  To  further  evaluate  the  versatility  of 
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the  method,  it  should  also  be  tested  using  other  types  of  structures 
such  as  different  aircraft  designs  and  components  or  building  structures. 
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APPENDIX  A 

STRAIN  GAGE  LOCATIONS  FOR  EXPERIMENTAL  PROGRAM 
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Three  wings  were  tested  in  the  experimental  portion  of  the  study. 
Surface  strains  were  measured  using  strain  gages  manufactured  by  Micro- 
Measurements  of  Romulus,  Michigan.  Two  types  of  gages  were  used:  EA- 
13-125AD-120  uniaxial  gages,  and  EA- 1 3-250RA- 1 20  three-gage  rectangular 
rosettes.  Figure  21  details  strain  gage  locations  for  all  three  tests. 
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The  rationale  for  sizing  rod  elements  representing  spar  caps  or  rib 
caps  was  presented  by  Jordan  (32,  33)-  All  symbols  in  the  following  sum¬ 
mary  of  his  presentation  refer  to  Figures  23  and  2k  which  were  extracted 
from  Reference  (32). 

Property  relationships  were  first  determined  for  a  structural  mem¬ 
ber's  cross  section.  The  distance  from  the  top  surface  of  the  section  to 
the  centroidal  axis  was  identified  as  ht<  Similarly,  the  distance  from 
the  centroidal  axis  to  the  bottom  surface  of  the  section  was  labeled  h. . 

D 

Maximum  bending  stresses  at  top  and  bottom  surfaces,  respectively,  were 


a 


T 


and 


(A. la) 


(A. lb) 


where  M  was  the  applied  bending  moment,  and  I  was  the  section  moment  of 
inertia  about  the  centroidal  axis. 

In  the  model,  the  rod  elements  were  assumed  to  be  point  areas  and 
were  positioned  at  the  top  and  bottom  surfaces  of  a  cross  section.  The 
rod  areas  were  sized  to  maintain  the  location  of  the  centroidal  axis  and 
the  value  of  I  for  the  actual  section.  Such  a  relationship  yielded 

\\  ■  Vb  (A-2) 

with  A_  and  A-  being  the  areas  of  the  top  and  bottom  rods.  The  bending 
T  D 

moment,  M,  in  the  model  then  became 

M  -  oTATht  +  oBABhb  (A. 3) 

Next,  using  actual  section  properties  and  desired  model  properties,  the 


relationships  of  Equations  (A. I),  (A. 2),  and  (A. 3)  were  combined  to  pro¬ 
duce 

AT  “  h  fh'  +h.  )  (A. 4a) 

t  t  D 

and 

AB  "  h.  (h  +h.  )  (A. 4b) 

D  t  D 

which  were  the  rod  element  areas. 

Skin  stiffeners  also  were  represented  by  rod  elements.  Initially, 
each  of  those  rods  had  the  same  cross  sectional  area  and  position  in  the 
structure  as  the  stiffener  it  represented.  However,  Models  A  and  B  used 
shear  panel  elements  to  represent  aircraft  skin,  so  the  membrane  capa¬ 
city  of  the  skin  panels  was  lost  because  shear  panel  elements  do  not 
represent  membrane  behavior.  The  membrane  capacity  of  the  skin  was  re¬ 
stored  to  models  A  and  B  by  adding  the  cross  sectional  area  of  each  skin 
panel  to  its  bordering  rod  elements.  This  is  illustrated  for  a  typical 
cross  section  in  Figure  24.  Such  modification  of  stiffener  rod  element 
areas  was  not  necessary  for  Models  C  and  D  because  they  used  membrane 
elements  to  represent  the  aircraft  skin. 
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(  A)  TYPICAL  CROSS  SECTION 


ROD  AREA  -  A,. 


SHEAR  PANEL' 


ROD  AREA  ■  Ag  '  1 

CROSS  SECTION 

MOMENT  *»  STRESS 
CONDITION 

(  B )  REPLACEMENT  SYSTEM 

Figure  23. 
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Check  the  first  element  in  the  model 
to  see  if  it  borders  the  failed  element 


Check  the  next  element  in  the  model 
to  see  if  it  borders  the  failed  element 


Check  margins  of  safety  for 
elements  in  the  next  group 


/  Did  \ 
/this  NASTRAN 
analysi s  use  the 


No 


5B 


L_J 

Make  no  changes  to  the 

■ 

case  control  deck 

Decrease  load.  Change  case  control 
deck  cards  for  LOAD,  SPC, 

MPC ,  and  LABEL 


Increase  load.  Change  case  control 
deck  cards  for  LOAD,  SPC, 

MPC,  and  LABEL 


List  failed 
elements 


List  elements  added  to 
the  damage  borders 


List  case  control  deck  cards 
for  LOAD,  SPC,  MPC,  LABEL  for 
next  NASTRAN  analysis 


Figure  25.  (Continued) 
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Figure  25.  (Continued) 
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APPENDIX  I 

PLOTS  OF  ANALYTICAL  AND  EXPERIMENTAL  DATA 
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Figure  26.  (Continued) 
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Figure  27.  (Continued) 
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(c)  Rod  Elemenl  564 

Figure  27-  (Continued) 
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(d)  Red  Element  576 

Figure  27.  (Continued) 
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Figure  28.  Comparison  of  Strains 
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Figure  28.  (Continued) 
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(c)  Rod  Element  564 

Figure  29 •  (Continued) 
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(a)  Load  Point,  Model  A  (without  torsional  stiffness  rods) 
Figure  31.  Single-Point  Displacements  for  Test  2C 
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(d)  Nod*  46.  Model  C  (without  torsional  stiffness  rods) 


Figure  31.  (Continued) 
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(a)  Load  Point.  Model  A  (without  torsional  stiffness  rods) 

Figure  32.  Single-Point  Displacements  for  Test  3B 
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(b)  Nod*  46.  Model  A  (without  torsional  stiffness  rods) 
Figure  32.  (Continued) 
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APPENDIX  J 

SUMMARY  OF  ANALYTICAL  RESULTS 
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TABLE  V 

SUMMARY  OF  RESULTS  FOR  TEST  1 , 
MODEL  0,  SIMPLE  DAMAGE 


Iteration 

No.  Load  Buckled  Elements  Overstressed  Elements  Failed  Elements 


1 

10 

68  118 

— 

— 

2 

12 

56  112 

— 

— 

3 

1A 

136 

— 

— 

A 

16 

SA  80  106  12A 

135  255 

529 

529 

5 

1A 

— 

— 

— 

6 

16 

— 

— 

— 

7 

18 

82 

— 

— 

8 

20 

111  117  132  I3A 
258 

— 

— 

9 

22 

6A  116 

117 

117 

10 

20 

— 

137  573 

573 

11 

13 

— 

119 

119 

12 

16 

— 

135  571 

135 

13 

1A 

— 

531 

531 

1A 

12 

— 

— 

— 

15 

1A 

— 

— 

— 

16 

16 

105  257 

— 

— 

17 

18 

85 

105  121 

105 

18 

16 

123 

533 

533 

19 

1A 

— 

— 

— 

20 

16 

— 

123 

123 

21 

1A 

— 

— 

— 

22 

16 

— 

— 

— 

23 

18 

— 

579 

579 

2  A 

16 

— 

— 

— 

25 

18 

— 

___ 

— 

26 

20 

--- 

— 

— 

27 

22 

52  62  89  107  125 

75  2A8 

2A8 

28 

20 

93  1 1A  1  A3  256 

107  125 

213  257 

125 
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TABLE  V  (Continued) 


Iteration 

No. 

Load 

Buckled  Elements 

Overstressed  Elements 

Failed  Elements 

29 

18 

126 

107  256  257  58l 

257 

30 

16 

— 

581 

581 

31 

14 

— 

126 

126 

32 

12 

— 

582 

582 

33 

10 

— 

— 

— 

34 

12 

— 

— 

— 

35 

14 

— 

— 

— 

TABLE  VI 

SUMMARY  OF  RESULTS  FOR  TEST  1, 
MODEL  D,  DETAILED  DAMAGE 


Iteration 

No. 

Load 

Buckled  Elements 

Overstressed  Elements 

Failed  Elements 

1 

10 

68  118 

— 

— 

2 

12 

56  112  136 

— 

— 

3 

14 

— 

— 

— 

4 

16 

80  135 

— 

— 

5 

18 

54  106  124  132 

— 

— 

6 

20 

116  117  13^  254 
255  258 

121 

121 

7 

18 

— 

117  137 

137 

8 

16 

Ill 

119 

119 

9 

14 

256  257 

577 

577 

12 

— 

— 

— 

11 

14 

123 

— 

— 

12 

16 

82  105  126 

123  531  9529 

123  9529 

13 

14 

— 

579 

579 

14 

12 

— 

— 

— 

15 

14 

125 

— 

— 

16 

16 

85  143 

105  125 

125 

17 

14 

— 

257 

257 

18 

12 

— 

— 

— 

19 

14 

— 

581 

581 

20 

12 

— 

— 

— 

21 

14 

— 

— 

— 

22 

16 

114 

— 

— 

23 

18 

62  64  89 

— 

— 

24 

52  93 

68  256 

68  256 

25 

18 

50  130 

— 

— 

26 

20 

55 

— 

— 

27 

22 

— 

126 

126 

aA  aa  aAa  pAa 


i» 
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TABLE  VI  (Continued) 


Iteration 

No. 

Load 

Suckled  Elements 

Overstressed  Elements 

Failed  Elements 

30 

20 

— 

— 

— 

31 

22 

— 

75 

75 

32 

20 

67  73 

55  57  473 

473 

33 

18 

97 

57  73 

73 

34 

18 

— 

55  471 

55 

35 

14 

53 

427 

427 

186 


TABLE  Vi  I 

SUMMARY  OF  RESULTS  FOR  TEST  2C, 
MOOEL  A,  SIMPLE  DAMAGE 


Iteration 

No.  Load  Buckled  Elements  Overstressed  Elements  Failed  Elements 


1 

10 

111 

135 

— 

-- 

- 

2 

12 

80 

— 

-- 

- 

3 

14 

49  81  82  133  138 

— 

-- 

- 

255 

4 

16 

68  1 

129  134 

109 

109 

5 

14 

114 

131 

127 

127 

6 

12 

132 

— 

-- 

- 

7 

14 

130 

— 

— 

- 

8 

16 

56 

— 

-- 

- 

9 

18 

50  79  107  108 

137 

137 

10 

16 

— 

119 

575 

575 

11 

14 

113 

123  125  143 

119 

129 

133 

135 

139 

119 

129 

579 

144 

573 

577 

579 

12 

12 

142 

103 

121 

131 

133 

135 

131 

577 

581 

139 

531 

567 

577 

581 

605 

615 

13 

10 

85  93  97  105  112 

113 

121 

125 

133 

135 

133 

582 

115 

116  117  118 

139 

140 

213 

223 

224 

124 

126  141  254 

225 

226 

227 

228 

535 

256 

257 

536 

567 

569 

573 

580 

582 

589 

591 

607 

611 

620 

625 

14 

8 

106 

97 

108 

115 

117 

121 

135 

528 

123 

125 

126 

135 

139 

140 

141 

223 

224 

225 

226 

227 

228 

257 

528 

534 

535 

536 

569 

571 

573 

578 

580 

589 

591 

609 

611 

625 

627 

629 

631 

632 

15 

6 

73 

97 

104 

108 

110 

113 

580 

115 

117 

121 

123 

125 

126 

128 

139 

140 

141 

144 

212 

213 

223 

224 

225 

226 

227 

228 

257 

530 

534 

535 

536 

569 

578 

580 

591 

593 

601 

Iteration 

No. 

Load  Buckled  Elements  Overstressed  Elements 

Failed  Elements 

15 

618  625  627  628  629 

(Cont . ) 

630  631  632 

16 

4  —  97  104  111  115  117 

118  121  123  124  125 

126  128  139  140  141 

142  212  213  222  223 

224  225  226  227  228 

249  257  532  534  535 
536  552  555  567  569 
578  589  591  593  595 

601  616  623  625  627 

628  629  630  631  632 

111  552  578 
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TABLE  VIII 

SUMMARY  OF  RESULTS  FOR  TEST  2C, 
MODEL  D,  DETAILED  DAMAGE 


Iteration 


No. 

Load 

Buckled  Elements 

Overstressed  Elements 

Failed  Elements 

1 

10 

56  68  80  112  114 
118  126 

— 

— 

2 

12 

106  132 

— 

— 

3 

14 

62  82  111  135  256 

— 

— 

4 

16 

49  50  108  124 

134  255 

- — 

— 

5 

18 

116  117  133  258 

68 

68 

6 

16 

— 

— 

— 

7 

18 

— 

66 

66 

8 

16 

79  81  113  129 

460 

460 

9 

14 

— 

— 

— 

10 

16 

— 

462 

462 

11 

14 

— 

— 

— 

12 

16 

— 

--- 

— 

13 

18 

61  105  107 

70 

70 

14 

16 

54 

— 

— 

15 

18 

52  115 

— 

— 

16 

20 

— 

56 

56 

17 

18 

63 

430 

430 

18 

16 

— 

— 

— 

19 

18 

— 

58  432 

58 

20 

16 

64 

432 

432 

21 

14 

— 

— 

— 

22 

16 

— 

— 

— 

23 

18 

— 

62  255  434 

62  255 

24 

16 

254 

64  78  81  82  187  434 

81  434 

25 

14 

— 

64  78  79  82  105  254 

436 

78  79  254 

26 

12 

55 

64  77  82  105  436  477 
479 

64  479 

Iteration 


No. 

Load 

Buckled  Elements 

Overstressed  Elements 

27 

10 

85 

77  82  105  436  477  481 

28 

8 

— 

— 

29 

10 

131 

54  80  428  475  477 

30 

3 

73  126 

52  75  82  428  438 

31 

6 

— 

55  473 

32 

4 

— 

— 

33 

6 

— 

427  429 

34 

4 

— 

— 

Failed  Elements 
77  436  481 

54  80  475  477 
52  75  428  438 

55 


427  429 


427  429 
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TABLE  IX 

SUMMARY  OF  RESULTS  FOR  TEST  3B, 
MODEL  A,  SIMPLE  DAMAGE 


Iteration 


No. 

Load 

Buckled  Elements 

Overstressed  Elements 

Failed  Elements 

1 

10 

— 

510 

510 

2 

8 

— 

— 

— 

3 

10 

— 

— 

— 

4 

11 

49  114  134  136 

110 

110 

5 

10 

130  132 

128 

128 

6 

8 

— 

— 

7 

10 

135 

— 

— 

8 

11 

— 

— 

— 

9 

12 

— 

— 

— 

10 

13 

116  133 

66 

66 

11 

12 

56 

— 

— 

12 

13 

111 

460  462 

460  462 

13 

12 

52  55  82 

424  432 

424  432 

14 

10 

80  113 

58 

58 

15 

8 

— 

528 

528 

16 

6 

— 

— 

— 

17 

8 

81  254 

434 

434 

18 

6 

— 

436 

436 

19 

4 

— 

— 

— 

20 

6 

79  256 

— 

— 

21 

8 

61  129  197  213 

52  438  482 

52  438 

22 

6 

62 

— 

— 

23 

8 

63  107  115  257 

109  127  437  582 

127  437  582 

24 

6 

— 

109 

109 

25 

4 

--- 

— 

— 

26 

27 

28 

29 

6 

8 

6 

4 

51 

54 

78  111 

60 

532 

78  111 

60 

532 

78  111 
60 
532 


78  III 
60 
532 
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TABLE  X 

SUMMARY  OF  RESULTS  FOR  TEST  3B, 
MODEL  D,  DETAILED  DAMAGE 


Iteration 


No. 

Load 

Buckled  Elements 

Overstressed 

Elements 

Failed  Elements 

1 

10 

56  68  112  114 

118  136 

— 

— 

2 

12 

132  256 

— 

— 

3 

14 

62  134 

— 

— 

4 

16 

49  80  116  135 

510 

510 

5 

14 

— 

110 

110 

6 

12 

130 

112  128 

112 

7 

10 

— 

566 

566 

8 

8 

— 

— 

— 

9 

10 

— 

— 

— 

10 

12 

— 

— 

— 

11 

14 

— 

114 

114 

12 

12 

50  82 

568 

568 

13 

10 

— 

— 

— 

14 

12 

— 

— 

— 

15 

14 

124 

— 

— 

16 

16 

258 

— 

— 

17 

17 

67  117 

— 

— 

18 

18 

111  254 

— 

— 

19 

19 

— 

530 

530 

20 

18 

— 

528  532  534 

574 

528  532  534 

21 

17 

55  64  73  97  107 
113  126  129  143 
255  257 

118  120  256 
574 

552 

120  256 

22 

16 

115  125 

116  117  118 
536  552  574 

138  257 

592 

117  118  536 

23 

14 

51  63  85  123  131 
133  141 

56  58  79  81 
115  119  137 
430  476  499 
574 

85  107 

257  279 
535  573 

73  107  115  279 
476  573  574 

23 


! 
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TABLE  X  (Continued) 


Iteration 

No.  Load  Buckled  Elements  Overstressed  Elements  Failed  Elements 

24  12  105  55  56  57  58  69  71  75  71  113  119  535 

78  113  119  197  429  571 

430  467  471  473  535 
571 

25  10  —  55  58  75  78  79  80  105  105  111  467  473 

111  121  197  429  430  569 

467  471  473  478  569 
575 

58  75  77  78  80  109  58  109  471  533 

197  429  430  471  533  567 

567 


26 
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